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1 Introduction 



The scenario where our world is associated with a brane embedded in a higher dimensional 
space-time with non-factorizable geometry [1-1] has attracted a lot of interest with the aim 
of solving several open questions in modern physics (see [5-8] for reviews). Most of examined 
models were realized as time independent field configurations. However, there have been con- 
sidered several braneworld models that assumed time-dependent metrics and fields [9-12]. Here 
we consider non-stationary braneworld scenario recently proposed in [13-15]. The braneworld 
is generated by 5D standing gravitational waves coupled to a phantom-like scalar field in the 
bulk. 

In our model [13-15] the bulk standing wave is bounded by the brane and the static part of 
the gravitational potential, which increases at the extra infinity. So there are some similarities 
with the experiments of [16], where neutrons were trapped between a reflecting plane at the 
surface of the Earth and the 'wall' generated by the gravitational potential of the Earth. We 
have discrete frequencies of standing wave oscillations just as the neutrons had discrete energy 
levels in those experiments. 

Another feature of the model [13-15] is the anisotropy of 4D part of the metric. Anisotropic 
braneworlds [17-23] have been considered when addressing various issues like anisotropy dissi- 
pation during inflation [ ], localization of test particles [20,21], or braneworld isotropization 
with the aid of magnetic fields [22]. As a general feature it has been established that anisotropic 
metrics on the brane prevent the bulk from being static [22,23]. 

A key requirement for realizing the braneworld idea is that the various bulk fields be localized 
on the brane. For reasons of economy and avoidance of charge universality obstruction [24] one 
would like to have a universal gravitational trapping mechanism for all fields. However, there 
are difficulties to realize such mechanism with stationary exponentially warped space-times. In 
the existing (l+4)-dimensional models spin-0 and spin-2 fields can be localized on the brane 
with the decreasing warp factor [3, 4], spin-1/2 field can be localized with the increasing warp 
factor [ ], and spin-1 fields are not localized at all [26]. In the case of 6D models it was found 
that spin-0, spin-1 and spin-2 fields are localized on the brane with the decreasing warp factor 
and spin-1/2 fields again are localized with the increasing warp factor [ ]. There exist also 
6D models with non-exponential warp factors providing gravitational localization of all kinds 
of bulk fields on the brane [28-32], however, these models require introduction of unnatural 
gravitational sources. 

In this paper we show that standing wave braneworld model [13-15] can provide universal 
gravitational trapping of zero modes of all kinds of matter fields in the case of rapid oscilla- 
tions of standing waves in the bulk. To clarify the mechanism of localization let us remind 
that standing electromagnetic waves, so-called optical lattices, can provide trapping of various 
particles by scattering, dipole and quadruple forces [33-35]. It is known that the motion of 
test particles in the field of a gravitational wave is similar to the motion of charged particles in 
the field of an electromagnetic wave [36]. Thus standing gravitational waves could also provide 
confinement of matter via quadruple forces. Indeed, the equations of motion of the system of 
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spinless particles in the quadruple approximation has the form [37]: 

—j— = Jpp = j^ D ^ R s , (1) 
as o 

where p M is the total momentum of the matter field and J a ^ 5 is the quadruple moment of 
the stress-energy tensor for the matter field. The oscillating metric due to gravitational waves 
should induce a quadruple moment in the matter fields. If the induced quadruple moment is 
out of phase with the gravitational wave the system energy increases in comparison with the 
resonant case and the fields/particles will feel a quadruple force, F M , which ejects them out of 
the high curvature region, i.e. it would localize them at the nodes. 

We proceed as follows. In Section 2 we recall basic ingredients of the 5D standing wave 
braneworld model [13-15]. In Section 3 we impose boundary conditions to choice exact form of 
background metric for matter fields and also list the approximations we use. Then, in Sections 
4, 5, 6 and 7 we demonstrate existence of zero modes of spin-0, -1, -1/2 and -2 particles on the 
brane, respectively. Short conclusions can be found in Section 8. 



2 Background solution 

5D standing wave braneworld [13-15] is generated by gravity coupled to a non-self-interacting 
scalar phantom-like field (different brane models with phantom fields can be found in [38-43]), 
which depends on time and propagates in the bulk. The action of the model has the form: 



S = J d 5 x^g 



1 ,„ _ . s 1 



(R - 2A 5 ) - ~g MN d. 



167rG 5 v ' 2 



M<pom 



(2) 



where capital Latin indexes refer to 5D space-time and G5 and A 5 are 5D Newton and cosmo- 
logical constants respectively. 

To avoid the well-known problems with stability, which occur with ghost fields, we can 
associate, for example, the ghost-like field <fi with the geometrical scalar field in a 5D integrable 
Weyl model [44-48]. In the Weyl model a massless scalar, either an ordinary or ghost-like 
scalar, appears through the definition of the covariant derivative of the metric tensor: 

Da9mn = 9mn9a<P ■ (3) 

This is a generalization of the Riemannian case where the covariant derivative of the metric is 
zero. The gravitational action for the Weyl 5D integrable model can be written as: 



Sw = J d 5 x^/g 



1 (R-2A 5 )-(6-509 MN d M <j>d m 



16nG. 



(4) 



where £ is an arbitrary constant. For £ = 11/10 the action (4) exactly coincides with (2). So 
we can start with a 5D Weyl model and require that we have Riemann geometry on the brane 
by assuming that the geometrical scalar field <fi is independent of 4D spatial coordinates and 
vanishes on the brane. 



3 



The definition of <fi via (3) avoids the usual instability problems of ghost fields since the 
geometrical fields have specific couplings with matter fields. They alter only lengths of vectors 
after parallel transport due to the assumption (3), but have not dangerous couplings like (pipip 
and it is known that the Weyl model is stable for any value of £. 

Situation is similar to the case with other geometrical structures. For instance, introduction 
of negative cosmological constants, or non-tensor character of Christoffel symbols does not lead 
to instability of particle physics models constructed in tangent spaces. 

The Einstein equations for the action (2) read: 

Rmn — ~9mnR = SkGbTmn — ^9mn ■ (5) 

We use the metric ansatz: 

ds 2 = e 2a|r| (dt 2 - e u dx 2 - e u dy 2 - e- 2u dz 2 ) - dr 2 , (6) 

where a is a constant, which corresponds to brane width. The peculiarity of the model (6) is 
that the brane, located at r = 0, possesses anisotropic oscillations and sends a wave into the 
bulk (as in [50,51]), i.e. the brane is warped along the spatial coordinates through the factors 
~ e u ( t < r \ which depend on time t and the extra coordinate r. 

The phantom-like scalar field <p(t, r) obeys the Klein-Gordon equation on the background 
space-time (6), 

4= d M (V99 MN d N <P) = e" 2a|r| - <j>" - 4a 8gn(r)<f>' = , (7) 

where sgn(r) is the sign function, overdots and primes mean derivatives with respect to t and 
r respectively, and determinant for our ansatz (6) is equal to 

V9 = e 4a|r| • (8) 
We further rewrite the Einstein - scalar field equations in the form: 

2 

Rmn = -d M ad N a + ~9mnK , (9) 
where the gravitational constant has been absorbed in the definition of the scalar field: 



a = v^ttCs <p . (10) 

It turns out that the system (9) is consistent only when the fields a and u are related with each 
other in the form [13-15]: 

The anisotropy of the metric (6) forces us to define the energy-momentum tensor of the 
brane with different tensions along different directions: 



M 
r N 



8(r) diag[A t , A^, X y , X z , 0] , (A x = \ y ) . (12) 
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The tensions A m (m = t,x,y, z) in general depend only on time. Taking into account (12) the 
field equations (9) change as follows: 



Rmn = -d M crd N a + -gMN A5 + 8nG 5 r MN 
where the reduced energy-momentum tensor, 

1 

r MN — T MN — -^9MN T , 

corresponds to the matter content on the brane, and takes the form: 

1 



tmn = g^Cr) diag [ 2X t + 2e~ u X x + e 2u X z , e u X t + X x - e iu X z 
e u X t + X X - e 3u X z , e~ 2u X t - 2e~ 3u X x + 2A 2 , ] 
The non-zero components of the Ricci tensor for the metric (6) read: 



-V 2a|T V + 4a 2 + 2a5(r) 



(13) 



(14) 



(15) 



Rtt 


= e Mr\ _. 


Rxx 


- R - e 2 

llyy C 


Rzz 


= e 2a\r\-2u 




3 a 
= ~2 U ~ 


R rt 


3. , 
= — uu . 
2 



2a\r\+u 



-e- 2alrl u - 4a 2 - 2ae(r)u' - 2aS(r) - -u" 



(16) 



Using the fine tuning: 



A, = 6a 2 



(17) 



the system of Einstein equations (13) reduces to a single ordinary differential equation, 



-2a\r\ 



r ' u — u" — 4a sgn{r)u = . 



In addition, for the brane tensions X m we have [ 

3a 



A, 

A.r 

A 2 
X r 



AttG ; 



-8(r) 

1 

3 «(*,o) 



e -2u(t,0) 

8vrG 5 
, 



(19) 



(6a + [u'])8(r) 



where u (t, 0) and [u'\ denote the value of u(t,r) and the jump of its first derivative at r = 0, 
respectively. 



The standing wave solution to (18) can be constructed by implementing the ansatz: 



u(t,r) = Csm{ut)f{r) , (20) 
where C and uj are real constants. From (18) we get the equation for the radial function: 

/" + 4a sgn{r)f + wV 2a|r| / = . (21) 
The general solution to this equation has the following form: 



f(r) = e 



-2a r 



AJ 2 [^e-^)+BYJ^e 
a \ I V a 



-a r 



(22) 



where A, B are arbitrary constants and J 2 and Y 2 are second-order Bessel functions of the first 
and second kind, respectively. 

3 Boundary conditions and averages 

As pointed out in [13-15] the ghost-like field a(t, r), along with the metric oscillations u(t,r), 
must be unobservable on the brane. Taking into account the relations (11) and (20) we can 
accomplish this requirement by setting the boundary condition for the function (22): 

/(r)| r=0 = 0. (23) 

Since J 2 and Y 2 are oscillatory functions, in the case of increasing (a > 0)/decreasing (a < 
0) warp factor for some fixed values of the constants A and B the function f(r) can have 
finite/infinite number of zeros. Thus the above boundary condition can be written in the form 
which quantizes the oscillation frequency, u, of the standing wave in terms of the curvature 
scale a, i.e. 

rr = x n , (24) 

\a\ 

where X n is the n th zero of the function f(r). Correspondingly, the nodes of standing wave, 
the points where the functions a(t, r) and u(t,r) vanish, can be considered as 4D space-time 
'islands', where the matter particles are assumed to be bound. 

For simplicity in this paper we assume that the constant A in the solution (22) is zero and 
the oscillatory metric function (20) has the form: 



u (t, r) = B sin (out) e~ 2a ^Y 2 f^- e" a|r| ) . (25) 



Physical solutions described by Bessel functions usually contain the first kind functions, J n , 
because of their regularity at the origin. However, in our case we can use even Bessel function 
of the second kind, Y n , since the argument in (25), e~ a ' r 'w /\a\, is always positive and becomes 
zero only for the case of increasing warp factor (a > 0) at r — > ±00. 
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To show explicitly pure gravitational localization of matter fields on the brane by the metric 
(6) we consider the case with the increasing warp factor a > 0. The oscillatory function (25) 
is zero at the position of the brane, r = 0, due to fine tuning (24), where X n now is one of the 
zeros of Y%. In this paper we explore the case when the Bessel function Yi has a single zero, i.e. 
we assume that 

- = X x » 3.38 . (26) 
a 

In the equations of matter fields the oscillatory function (25) enters via some exponential 
functions: 

{buf 



E^> < 27 > 



n=0 



where b is a constant. We suppose that the frequency of standing waves u in the oscillatory 
metric function u(t,r) is much larger than frequencies associated with energies E of particles 
on the brane, i.e. 

u > E . (28) 

In this case we can perform time averaging of oscillating exponents in the equations of matter 
fields. 

Using the expression: 



(m = 2n + l) 



— / [sm(uJt)} m dt = < 2 -™m\ (m-2n) ( 29 ) 

2tt J [ [m - In) 

for the time averages of oscillating exponents (27) we get the simple formula [52]: 

<«-> = E^ = W. (30) 

where Iq is the modified Bessel function of the zero order of the argument: 

Z = \bB\e- 2aM Y 2 (^e" a|r| ) . (31) 
Let us also display some equalities for time averages of oscillatory functions [53]: 

<«> = M = ^W£«^=0. (32) 



t dt/ \dt 

where the prime denotes the derivative with respect to the extra coordinate r. 

Now we are ready to consider localization problem of various matter fields on the brane 
within the standing wave braneworld model with background metric (6), where the oscillatory 
function u(t,r) has the form (25). 
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4 Localization of scalar fields 

We start with the problem of localization of massless scalar fields defined by the 5D action [53]: 



S * = ~IJ V9dx 4 dr g MN d M ®d N $ . 



(33) 



Corresponding Klein-Gordon equation, 



9m {Vgg MN d N <5>) = , (34) 

can be written as: 

[d 2 - e~ u {dl + d 2 ) - e 2u d 2 z ] $ = e 2a|r| (e 4a|j W)' . (35) 
We look for the solution of this equation in the form: 

$ (t, x, y, z, r) = r) X (x, y)£(z) . (36) 
This separation of variables transforms (35) into the system of the equations, 

{d 2 x + d 2 ) X +{pl+p 2 y )x = 0, 

dit + vli = 0, (37) 
9^ + M + P 2 y )e- U + p 2 z e 2u ] * = e 2a ^ (e Aa ^')' . 

On the brane, where u ~ and anisotropics are absent, the parameters p x ,Py and p z can be 
regarded as momentum components along the brane. Their exact physical meaning is not clear 
in general case since (37) are equations with variable coefficients. 
Now we separate variables in the last equation of the system (37), 

tf(t,r) = e iE \{r) , (38) 

where the energy E for 4D massless scalar modes on the brane obeys the dispersion relation: 

E 2 =p 2 x +p 2 y +p 2 z . (39) 

As it was mentioned in previous section, when the frequency of standing waves u is much larger 
than frequencies associated with the energies of the particles on the brane, 

u > E , (40) 

we can perform time averaging of oscillatory functions. Then for the extra factor q{r) of the 
scalar field wave function we obtain the equation: 

(e 4a ! r V) / — e 2a ^K(r)q = , (41) 
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which contains the time-independent factor: 

K(r) = «e-> - 1) ( P l+P 2 y ) + ((e 2u )-l)p 



(42) 



It is more convenient to put (41) into the form of an analogue non-relativistic quantum 
mechanical problem by making the change: 



-2a\r 



'l-0 s (r) 



For ip s { r ) we fi n d: 
where the function 



<-f/ s (r)^ = 



(43) 
(44) 



U s (r) = Aa5{r) + 4a 2 + e _2a|r| i<:(r) (45) 

is the analog of non-relativistic potential. FIG. 1 shows behaviour of U s (r) in the case defined 
by (26). 




Figure 1: The effective bulk potential (45). 



To study general behaviour of extra part of the scalar field zero mode wave function we 
explore (44) in two limiting regions, far from and close to the brane. The function K(r), which 
describes oscillatory properties of standing waves, has the following asymptotical forms: 



*(r)U, 



~ r , 
~ const . 



(46) 
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So close to the brane, r — > 0, the equation (44) obtains the form: 



- [Aa5{r) + 4a 2 ] tp s = , (47) 

with the solutions 

iWUo~^ H • (48) 

At the infinity the equation (44) again reduces to (48) (but without delta- function of course) 
with the two solutions, 

^MUo-e^l . (49) 

Resemblance of the equations far from and close to the brane is not surprising since standing 
waves in the bulk are bounded by the brane at r = and by the gravitational potential at the 
infinity, i.e. these two regions correspond to the nodes of the waves. 

To have normalizable zero mode we impose the boundary conditions: 

^(r)U = 0, 

fWUoo = °- ( 5 °) 
Then using the definition (43) and the solutions (48) and (49) we find: 

<^( r )lr^o ~ const, 

*MU» ~ e" 4aM - (51) 

So q(r) has a maximum on the brane and falls off at the infinity as e _4a ' r '. In the action of 
scalar fields (33) the determinant (8) and the metric tensor with upper indices give the total 
exponential factor e 2a ' r ', which obviously increases for a > 0. This is the reason why in the 
original brane models [3, I] scalar field zero modes, with the constant r-depended extra part 
of the wave function, can be localized on the brane only in the case of decreasing warp factor 
(i.e. a < 0). In our model the extra part of wave function (51) is not a constant, moreover, for 
a > it contains the exponentially decreasing factor e _4a ' r L For such extra dimension factor 
the integral over r in the action (33) is convergent, i.e. 4D scalar fields are localized on the 
brane. 

5 Localization of vector fields 

For simplicity in this paper we investigate only the U(l) vector field [54], the generalization to 
the case of non-Abelian gauge fields is straightforward. The action of vector field is: 

S v = -\J d 5 x^ g MN g PR F MP F NR , (52) 

where 

F M p = d M A P - d P A M (53) 

is the 5D vector field tensor. 
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Incidentally, in 5D space-time implementation of the pure gravitational trapping mecha- 
nism of vector fields on the brane remains most problematic. In the original brane models 
extra dimension parts of the vector field zero modes are constant functions. So, taking into 
account that in the action integral (52) the extra coordinate parts of the determinant (8) and 
two metric tensors with upper indices cancel each other, the vector field zero modes are not 
localizable on the brane for any sign of a. Because of this fact there was proposed some non- 
gravitational trapping mechanisms, for example [55]. Here we want to show that the standing 
wave braneworld metric (6) provides pure gravitational localization of vector field zero modes 
on the brane. 

The action (52) gives the system of five equations: 

^pd M (Jg g MN 9 PR F NR ) = . (54) 

Let us seek for the solution to the system (54) in the form: 

A t (x c ) = v{r) a t (x u ) , 
A x (x c ) = e<^v{r) <x v ) , 

A y (x c ) = e u{t ' r) v(r) a y (x u ) , (55) 
A z (x c ) = e- 2u{t ' r) v(r) a z (x u ) , 
A r (x c ) = , 

where u(t, r) is the oscillatory metric function (25), a IM (x u ) denote the components of 4D vector 
potential (Greek letters are used for 4D indices) and scalar factor v(r) depends only on the 
extra coordinate r. The last expression in (55) is in fact the 5D gauge condition. 
We require existence of flat 4D vector waves localized on the brane, 

a„ {x u ) ~ £ ^ Et+p * x+p y y+p * z) , (56) 

where E, p x , p y , p z are components of energy-momentum along the brane. Solutions of this kind 
exist only on the brane, where u = 0, and in the case when u is much larger than frequencies 
associated with the energies of the particles on the brane E. For such high frequencies of bulk 
standing waves we can perform time averaging of oscillatory functions in (54). 

Taking into account the equalities (32) time averaging of the fifth equation of the system 
(54), 

d a {fPA'p) = , (57) 

gives the Lorenz-like gauge condition: 

g afi d a A fi = v^d^ = , (58) 

where 7] a p denotes the metric of 4D Minkowski space-time. The equation (58) together with 
the last expression of (55) can be considered as the full set of imposed gauge conditions. 

Performing time averaging and making use of (32) and (58), remaining four equations of 
the system (54), 

d, (g^g Pa F Sa ) --^{Vg 9 Pa K)' = , (59) 
\J9 
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reduce to: 

v g aS d a d & a p + e~ 2a|r| (e 2a|7 V)' a p = . (60) 
In the case of the zero mode (56), the system (60) gives the single equation for v(r): 

(e 2a ' r 't/)' — K(r)v = , (61) 
where K{r) is defined by (42). By making the change: 

v (r) = e- a ^ v (r) , (62) 

as in the case of scalar field (44), we put the equation (61) into the form of a non-relativistic 
quantum mechanical problem: 

<(r) - £/ v (r)W(r) = , (63) 

where the potential 

U v (r) = 2a5(r) + a 2 + e - 2a|r| if (r) (64) 

differs from the analogues potential for scalar fields (45) (see FIG. 1) only in constant factors 
in the first and second terms. So the solutions to (63) close to and far from the brane are: 

^(r)U» ~ e ±a H. (65) 



Taking into account the definition (62) and imposing boundary conditions analogues to (50), 
for the extra dimension factor of the vector field zero mode wavefunction we get: 

^( r )lr^o ~ const, 

^(OUoo ~ e-**"*-!. (66) 

It's easy to see that the extra factor v(r) of the vector field zero mode wave function has 
maximum on the brane and falls off at the infinity as e _2a ' r L For such zero mode solution 
integrals over r in the action of vector fields (52) are convergent, therefore 4D vector fields are 
localized on the brane. 

Indeed, using the ansatz (55) and the equalities (32), time-averaged components of the 
vector field tensor (53) can be written in the form: 

(F; (x c )) = {g vP F rP ) = e~ 2a ^v' a v (x a ) , 

(F?{x c )) = (g^F^) = e- 2a \%[fU(x a ) + MH((e bu -l)dea a )] , (67) 

where 

fiw (x a ) = d^a u - d u a^ (68) 

is the 4D vector field tensor and MfJ represents the sum of terms of the type (e bu — l) d u a^. 
Since the functions v 1 and (e feu — l) in (67) vanish on the brane, the 5D vector lagrangian at 
r = gets the standard 4D form: 



= - V(0) fit, ■ (69) 

r=0 4 
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Of course in general the 5D lagrangian of vector fields is more complicate. For instance, it 
contains the mass term (v'F^) 2 ~ {v') 2 a a a a , i.e. the zero mode vector particles are massless 
only on the brane and acquire masses in the bulk. This fact can be considered as the alternative 
mechanism of localization. 



6 Localization of massless fermions 

In this section we investigate the localization problem for massless fermions within the standing 
braneworld model [56]. For Minkowskian 4x4 gamma matrices {{'J a ,'j 13 } = 2rj a ^) we use the 
Weyl basis, 



I \ , / -a 



( /()r \ a 1 



>y" = yyt l 'y z 



I 
-I 



where I and a 1 (i = x, y, z) denote the standard 2x2 unit and Pauli matrices respectively. 

Let us recall that four-component columns represent fermions in 5D, and that 5D gamma 
matrices can be chosen as: 

Y A = h A A Y A , {T A ,T B } =2g AB , 

r^= (y, 7 *,7^, 7 V 7 5 ) , (70) 

where A, B, ... refer to 5D local Lorentz (tangent) frame. So according to (70) the curved-space 
gamma matrices T A are related to Minkowskian ones as: 

7* , 

-«/2 7 y ? ( 71 ) 



The fiinfbein for our metric (6), 

h\ = (e a ' r ', e a ' r ' +n/ ' 2 , e a ' r ' +n//2 , e a ' r '~", l) 
h AA = g AB h A , hi = VAS h BA , h AA = VAB hl (72) 

is introduced through the conventional definition: 

9ab = VABhihs • (73) 

The 5D Dirac action for free massless fermions can be written as: 

S = J d 5 x^ (x A ) T M D M q (x A ) , (74) 
13 









= e~ a 


T y 


= e~ a 


V z 


= e~ a 


r r 


= i'j 5 



where covariant derivatives are defined as follows: 



D A = d A + -nfTzTc . (75) 



In the last expression VL^ N denotes the spin- connect ion: 

1 



C)MN _ _Cl NM 



h NM (d M h% - d N hl) - h NR (d M h% - d N h$) ~ (76) 
-h A M h pM h? N (d P h QA -d Q h pA ) 



The non-vanishing components of the spin-connection in the background (6) are: 

flf = - (e a|r| ~ n )' , (77) 

- d (e u/2 ) 
o xt — Of* — v ' 

Ux ~ u y ~ dt ' 



n 



zt 



d(e~ u ) 



dt 

The corresponding to (74) 5D Dirac equation reads: 

iT A D A ty = i (VD^ + T r D r ) * = . (78) 

For the bulk fermion field wave function we use the chiral decomposition: 

* (x\ r) = ^ L {x v ) A(r) + fa {x v ) p(r) , (79) 

where A(r) and p(r) are extra dimension factors of the left and right fermion wave functions 
respectively. We assume that 4D left and right Dirac spinors, 

= -i>L , l^R = , (80) 

correspond to zero mode wave functions, i.e. they satisfy free Dirac equations: 

i-fd^L = i-fd^R = . (81) 

Apart from the massless states ipL and i/jr, the 5D Dirac equation also have solutions cor- 
responding to massive fermions. In the single brane models the masses of the bounded massive 
states are typically of order of the energy scale a, characterizing the brane as a topological 
defect in higher-dimensional space-time. These states are very heavy and we do not consider 
them here. 
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The solutions of (81) in our representation can be written in the form: 

ipi, (x v ) = ( ? J e-^-P^-Pyy-P^) } (82) 



x L / 

where the constant 2-spinors L and R satisfy the relations: 

(E + a%) L=(E- a l Pi ) R = . (83) 

When the frequency u of standing waves in the background metric (6) is much larger than 
frequencies associated with the energies E of the fermions on the brane we can perform time 
averaging of the oscillatory functions in the Dirac equation (78). Time averages of the Dirac 
operators are: 



(<r*A> 


= ie~ alrl 


7*9t - \, sgn(r)^ 5 , 


iiY x D x ) 


= ie~ alrl 




(iT y D y ) 


= ie~ alrl 




(iT*D z ) 


= ie~ alrl 


(e u ) i z d z - l -a sgn(r)j 5 , 


(iT r D r ) 







34) 



and the equation (78) can be written in the form: 

i [7% + {e u ' 2 ) (^d x + i«d y ) + (e- u ) 7 Z «9 Z ] * 



= e a|r| 7 5 [2a sgn(r) + d r ] ¥ . (85) 

Using the solutions of free equations (82) and the relations (83), the equation (85) can be 
rewritten as the system: 



-e a|r| [2a sgn{r) + d r ] <rVP*(r) \ ( p{r)R 

-(jiPiir) e Q l r l [2a sgn{r) + d r ] J \ \{r)L 

Here we have introduced V -dependent momentum' Pi(r): 



. (86) 



P x {r) = ((e- u / 2 )-l)p x = [I (Z(r)/2)-l]p x , 

Py(r) = {(e^ 2 )-l)p y = [I (Z(r)/2)-l}p yj 

Pz(r) = ((e u )-l) Pz = [I (Z(r))-l]p z , (87) 

P\r) = P 2 + P 2 + P 2 . 



where Z{r) is defined by (31). 
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From the second equation of the system (86) it is straightforward to find 



p ( r )R = e a l r 'l^^ [2a sgn{r) + d r ] X(r)L . (88) 
P A (r) 

Inserting (88) into the first equation of (86) and multiplying the result by a 1 Pi, we receive the 
second order differential equation for the function A(r): 



A" + 



5a sgn(r) — 



P' 



A' + 



P' 



Aa5{r) + 6a 2 - 2a sgn{r) — - P 2 e " 2a|r| 



P 



A = 



S9) 



Now, as in the cases of scalar and vector fields, let us investigate this equation far from and 
close to the brane. 

Close to the brane, r — > ±0, the V -dependent momentum' (87) behaves as: 



P(r) 



Ar 2 + 0(r 3 ) , 



where A is constant, and the equation (89) takes the following asymptotic form: 



A" + 



5a sgn(r) 



A' + 



4a<5(r) + 6a 2 sgn(r 



A = . 



This equation has the unique nontrivial solution: 

A(r)U ±0 = Ce~ 2a ' r ' 



(90) 



(91) 



(92) 



where C is a constant. 

As it follows from (88) and (92), in our setup the right fermionic modes are absent on the 
brane: 

p(OU±o = 0. (93) 

Such different behavior of the left and right massless fermions on the brane is not surprising, 
since in our model the effective mass term in (85) is of 7 5 -type, 



m(r) = 2a , y 5 sgn(r)e 



a\r\ 



with the gap: 



\m(r) — m(—r)\ = Aa^ e 



5 „a\r\ 



(94) 



(95) 



In the second limiting region r — > ±oo the function P'/P vanishes and the equation (89) 
obtains the asymptotic form: 



A" + 5a s^(r)A' + 6a 2 A = 



with the solution: 



K r )\r^±oo ~ e 



-3a I r 



(96) 
(97) 
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Using (97) from the relation (88) we find the asymptotic behavior of the extra factor of the 
right fermion wave function: 

p(r)\ r ^ ±00 ~ e -^\r\ . (98) 

So in our model the extra dimension part A(r) of the bulk left spinor wave function (79) 
has the maximum at the origin, 

X(r)\ r=0 = C, (99) 

falls off from the brane, and turns into the asymptotic form (97) at the infinity. When r — > oo 
the determinant (8) in the action integral for 5D fermions (74) increases as e 4a ' r ' for a > 0. 
However, extra dimension factor of left fermions, according to (97), contribute e _6a ' r ' and overall 
r-depended part in (74) decreases as e _2a ' r '. So in the case of left fermions the integral over r 
in (74) is convergent, i.e. zero modes of left fermions are localized on the brane. 

At the same time, according to (93) and (98), right fermion zero modes does not exist on 
the brane and the extra dimension part p(r) of their wave functions at the infinity decreases as 
e -2a|r|_ therefore for right fermions integral over r in (74) diverges and zero mode wavefunctions 
of right fermions actually are not normalizable. 

7 Localization of gravitons 

Finally, let us deal with spin— 2 graviton. We shall consider the metric fluctuations: 

ds 2 = e 2a|r| (g^ + V) dx»dx v - dr 2 , (100) 

where g^ u is the metric tensor of the 4D part of (6): 

g, u = (l,-e u ,-e u ,-e- 2u ) . (101) 

Let us suppose that the frequency of standing waves in the bulk u is much larger than 
frequencies of gravitational waves on the brane, i.e. we can replace the oscillatory functions 
e u{t,r) j n (ioi) D y their time averages (30). So, since the brane is placed in a node of the standing 
wave, close to r = we can use the approximation: 

( e «(*.»0) + . (102) 

This means that the functions (u) can be regarded as r-depended additive terms of h^ v in flat 
Minkowski space. Then the equations of motion for the fluctuations h^ u are found to be 

4= q m (Vgg MN dNh^) = o . (103) 

It turns out that these equations are equivalent to the equation of motion of a spin-0 scalar field 
(34) if we replace $ with h^ u . Accordingly, the localization problems for spin-2 graviton field 
and spin-0 scalar field, considered in the Section 4, are similar, and gravitons are also localized 
on the brane. 
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8 Conclusion 



In this letter we have investigated the localization problem of all kinds of fields (scalar, vector, 
spinor and tensor) within non-stationary braneworld scenario [13-15], where the braneworld 
is generated by 5D standing gravitational waves coupled to a phantom-like scalar field in the 
bulk. It is noteworthy that the trapping of fields on the brane in our model has universal and 
purely gravitational nature. The conspicuous feature of our model, as opposed to earlier static 
approaches with decreasing warp factors [ , ], is the use of the metric ansatz with increasing 
warp factor. By explicit calculations we have shown that scalar, vector and tensor field zero 
modes are localized on the brane. In the case of fermion fields we have found that the left 
fermion zero mode is localized on the brane, while the right fermion zero mode does not exist 
on the brane and its wavefunction actually is not normalizable. 

In our opinion one can use the restriction (26), which controls number of nodes of the bulk 
standing waves in our model, i.e. number of parallel 4D space-time islands, to address the 
old problem - the nature of flavor. In this regard our scenario has to be studied in detail 
using numerical solutions to the equations of matter fields in the case of interacting fields. 
Corresponding investigations are currently in progress. 
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